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Abstract 

Ionic-pair equal-time displacement correlations in relaxor lead magnesium 
niobate, Pb(Mg 1 / 3 Nb 2 /3)C'3, have been investigated at 300, 270, and 240K 
in terms of an x-ray diffuse scattering technique. Functions of the distinct 
correlations have been determined quantitatively. The results show the sig- 
nificantly strong rhombohedral-polar correlations regarding Pb-O, Mg/Nb-O, 
and O-O' pairs. Their spatial distribution at 300K forms an ellipse or a sphere 
with the diameters of 30-804, and the sizes are reduced to 30-404 on cooling 
through local condensation of thermal lattice fluctuations. This observation 
of local structure in the system proves presence of the polar microregions in 

the paraelectric state which leads to the dielectric dispersion. 
PACS numbers: 77.80.Bh, 61.10.Eq, 64.70.Kb, 77.84.Dy 
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Pb(Mgi/ 3 Nb2/3)0 3 (called PMN) is a representative system of relaxors, of a special kind 
of ferroelectrics exhibiting a diffuse cubic-to-rhombohedral phase transition and dielectric 
dispersion [1,2] near 270K. The crystal structure, based upon a simple cubic perovskite often 
represented by ABO3, includes two different ionic species at site B, Mg 2+ and Nb 5+ , which 
are known to make 1:1 chemical short-range order (SRO) along the <111> [1,2]. This local 
atomic nature has been suggested to cause the unusual phenomena, but a microscopic origin 
of the mechanism is yet to be fully understood. 

The ordering nature of PMN is expected to cause two kinds of lattice distortion. The 
one kind is of an electrostatic origin. B-site occupation of the two ions gives a bulk-averaged 
valence of +4, and this value electrically neutralizes the whole system with the other species, 
Pb 2+ , 0\~ , C>2~ ) an d Olr (Subscripts of three oxygens mean the cell faces where each species 
is centered. The faces are perpendicular to x,y, and z,respectively). However, because of 
the SRO, the valence averaged in the ranges becomes +3.5 and hence ~5+ in the Nb-rich 
surroundings. This charge misfit is suggested to produce lattice distortion electrostatically 
[3]. The other kind lies in a mechanical origin, a difference in ionic sizes between the B- 
site ions [4,5], called size-effect, which also leads to the same consequence as that of the 
former kind. Both effects can cause local destruction of the cubic symmetry giving the 
permanent dipoles whose presence enables us to evaluate models previously proposed to 
explain a mechanism of the relaxor behavior [6,7]. Therefore, it is very important to directly 
observe local ionic displacements in PMN to obtain experimental evidence of existence of 
the polar correlations, which motivated us to carry out a measurement of x-ray diffuse 
scattering. There are preceding experiments performed based upon the same purpose by 
means of NMR [8] to investigate the polar time correlations, and by a transmission electron 
microscopy (TEM) [9] and inelastic neutron scattering (on PZNT which is the different 
system) [10] to observe the spatial correlation. Our x-ray diffuse scattering work will be 
compensative for these previous studies, giving the lengths and the precise symmetry. 

An x-ray (also neutron and electron) intensity scattered by atoms in a crystal, defined as 
I(qi) where q\ is a scattering vector, consists of the components given by an averaged crystal 
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structure and local atomic arrangements, lBragg(qi) and lD%ff(qi), respectively. The former 
gives a Bragg intensity, and the latter a diffuse scattering intensity which is tremendously 
weaker than the other component. In the current case, only Ioiffiqi) is dealt with, and is 
further divided into the contributions from the local displacements (including phonons) and 
the SRO, Idisp{Qi) and Isro(Qi)- F° r mathematical expression of the two intensities, the 
simple way is to adopt sums of plane waves [11] which represent the atomic configurational 
disturbance. For example, the displacement of the species 7 at the s th unit cell, u sl , is 
written in the following form: 

« S 7 = EE u jkA ex P\ 1 ^ ■ {Rs + r 7 )] , (1) 
k » 

where i is an index of real-space Cartesian coordinates, u ^ is the displacement-wave am- 
plitude of wavevector k allowed in the first Brillouin zone, is a unit vector, and -R s +r 7 
is the atomic position. The concentration is also given in the same way, and has the form 
similar to Equation (1). If these plane- wave sums are used for the expression, the intensity 
equations obtain the simple forms in electron units (e.u.). We show only IoispiQi) as below: 

iDisp(qi) = Jv2 EEEE A/y exp[i(k + gl) • (r 7 - ry )] 

7 7' i V 

(ei ■ qi) {e v ■ q x ) < u^y > 8(k + q x - 2nH hkl ) , (2) 

where N is the number of unit cells considered, / 7 is the x-ray atomic scattering factor 
[12] including the temperature factor [13], and 2nH hk i is a Bragg reciprocal lattice vector 
(remove (e^ • qi)^ ■ q\) term to gain the form of IsRo(qi))- A bracket of amplitude product 
< u jki u yki' > originates from a fact that the correlations which we observe are spatially 
averaged, < • • • >. The products are only unknowns in Eq.(2) and their values are unique 
for all vectors q\ with unique k. For PMN, there are 17 unknowns including the SRO terms 
(The number seems small, but is a result from considerable degeneracy among the terms 
because of space group Pm3m.). Therefore, if at least 17 total scattering intensities observed 
for universal k are obtained, the unknowns can be determined through linear algebra. The 
computed values for all allowed k provide the ionic-pair equal-time displacement correlation 
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functions according to the following definition: 

< u sji uwv > = J2 < u -yki u Yki> > ex P[ l k ■ ((Rs + r-r) - (0 + ry ))] , (3) 
k 

where R s > is set at the real space origin, s'=0. A plot of Eq.(3) as a function of R s represents 
strength and distribution of the 7-7' correlation required for this study. 

An x-ray diffuse scattering measurement was performed at room temperature (300K) in 
the air and below 300K with a DISPLEX cryostat on self-flux-grown single crystal PMN 
whose size is ~6x4x2mm 3 . A surface normal of the largest face is almost parallel to 
the [001], and with respect to this face, the bisecting-26* diffraction was done on a 4-circle 
diffractometer (Hiiber) with a scintillation counter. In this way, x-ray absorption by the 
sample is independent of q\. The conventional Cuka rotating-anode generator (Mac Science) 
was used and operated at 15kW. X-rays were monochromated by a pyrolytic graphite (002). 
A preliminary measurement of the Bragg intensity distribution showed that the volume- 
distribution extent is ~(0.01) 3 in reciprocal lattice units (r.l.u. 3 ). Therefore, minimum 
I k I was determined as 0.02 r.l.u. to avoid the Bragg contamination. The intensities 
were collected in reciprocal lattice volume with changing qi by increment 0.02 r.l.u. (i.e., 
N=50 3 ). The volume contains about 5000 independent q\s. The points consist of about 160 
qi sets, and each set includes at least 17 points. The data acquisition time at each point 
was 5 minutes, and the complete measurement took about 3 weeks. A factor to convert the 
observed intensities to those in e.u. was determined from the integrated Bragg intensities 
of an aluminum standard powder sample [14]. Experimental conditions and means for this 
measurement are the same as those for the diffuse scattering measurement. The contribution 
of two kinds of unwanted scattering, air scattering and Compton scattering, was subtracted 
from the observed intensities through the measurement and the calculation [15], respectively. 

Figure 1 (a) is a contour map of the observed e.u. intensities at 300K in the 0th (100) 
reciprocal lattice plane. A portion of the collected data was used for this mapping. Ridges of 
the intensities are seen around all Brillouin zone centers along the [Oil]. The tiny ridges along 
the [011] also appear around some zones, indicating <fi dependence on the distribution. This 
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feature essentially agrees with that observed in previous x-ray work [3] which suggested the 
origin from the lattice distortion making competing interactions with the <011> transverse 
optic phonons. The current results show the significant intensities around the zone centers 
whose three digits of the indices h, k, and 1 give the following relationship, h+k+l= even 
integer. The point is that these scattering intensities are nearly isotropic close to the zone 
centers. Such zones were expected to exhibit a little to no ANISOTROPIC diffuse scattering 
in the previous work [3], which is consistent with the present results. Emergence of these 
nearly isotropic intensities suggests that the optic mode which we are dealing with is not 
purely optic [16]. The impure component may arise from uncollaborative interactions of the 
charge misfit with the pure mode. The observed data were analyzed based upon Eq.(2), 
and < u ^ u > terms were determined quantitatively. The intensities were synthesized 
using the determined parameters. The recalculated map is represented in Fig. 1 (b), and 
the figure reproduces the observed data very well giving the R factor ~2%. 

Results of the parameter values indicate that the three distinct displacement amplitude 
products, < u ph%x u Q1%y >, < u B%x u Ql%y >, and < u ol%y u Q2%z >, are significant only around 
| A; |=0 (species B means the B-site ion, Mg 2+ / Nb 5+ .). Their values are largely negative, 
positive, and positive quite near the zone center, and oscillate almost in phase with each 
other with increasing | k |. Since their absolute values are 20-10 5 times bigger than all 
the others, only these three terms are required to describe the local ionic structure. Its 
qualitative representation is given below. 

The negative values of < u Pb ^ x u Q1 ^ y > means a product of the shifts of Pb and 1 
along, e.g., -x and +y, respectively. The shifts are illustrated on the left-hand side in 
Figure 2 (a). Because of the simple-cubic-perovskite symmetry, < u pb j: x u 01 j: y > is equiv- 
alent to < u Pb j: x Uoikz > - Involving 2 and 3 , the further symmetry consideration in- 
dicate more equivalencies such as < u Pb ^u ol ^ y >=< u pbiy u ^ x >=< u p ^ y u 02 ^ z >=< 
u Pbkz u 03kx > =< u pbkz u 03ky > - Therefore, superposition of all the equivalencies gives the 
resultant shifts along the [111] for Pb, the [Oil] for O u the [101] for 2 , and the [110] for 3 . 
This picture is seen in the otther two figures of Fig. 2 (a. Similary, as shown in the right- 
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hand figure of Fig. 2 (b), positive < u B ^ x u ol ^ y > implies a product of the shifts of B and 
0\ along +x and +y, respectively. The equivalencies can be written as < u B k x u oik y > =< 

U Bkx U 01kz > = < U Bky U 02kx > = < U Bky U 02kz > = < U Bkz U 03kx > = < U Bkz U OSky >■ A reSult 

of the superposition provides the shift of B along the [111] given in the other figure of Fig. 2. 
A set of the oxygen shifts shows no contradiction with the outcome in the Pb-0 case. The 
rest, positive < Uoiky u 02kz > ' ^ s displayed in Fig. 2 (c), and the equivalencies are given as 

< U 0\ky U 02kz > = < u O\kz U 02kx > = < U 0\ky U 03kx > = < U 0lkz U 03ky > = < U 02kx U 03ky > = < 

u 02kz u 03kx > - The resu hant schematic shifts are exactly the same as those of the two cases 
above. In addition, the parameter value for the concentration is also significant only around 
£=1/2<111> with which oscillation of the concentration, i.e., the 1:1 SRO, is exhibited. 
But we are not concerned with the SRO in this paper since the concentration correlation 
has been pretty much studied so far [9,17]. 

A combination of all the shifts derived from the three main correlations above leads to 
the local structure shown in Fig. 2 (d). Two pyramids confronting each other whose vertices 
are occupied by Pb are both translated along the [111] from the ideal cubic positions given at 
corners of the shaded box, and maintain point group 3 which the cubic structure also has. On 
the other hand, the other two B-0 pyramids between the two of Pb linking at B realizes the 
lower symmetry, 3, because of the displacements of B and O. Therefore, this local structure 
totally forms rhombohedral polar symmetry which was suggested by previous x-ray and 
neutron 'averaged' crystal structure analysis at low temperature [13]. The similar analysis 
was also performed by neutron diffuse scattering at single k [18] near the zone centers, and 
the results show general agreement with our local structure. To our knowledge, existence of 
the polar rhombohedral correlation has crystallographically been proved directly from the 
diffuse scattering analysis for the first time, and there is no contradiction with results of the 
NMR work [8]. 

Quantitative representation of the correlations was done by computing < u s pb x uooiy >, 

< u S BxUooiy >, and < u sOly u o2z >, according to Eq. (3). They are given in Figures 3 
which shows the plots with respect to intercell distance, | R s |, along (a) the [111], (b) [110], 
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and (c) [100] directions. For all figures, values of < u s oi y uoo2z > are substantially ten times 
larger than those of the other two functions. This difference is suggested to originate from 
phonons which reveal larger displacement amplitudes for lighter elements. Since our x-ray 
scattering does not distinguish phonons from static displacements, the functions also involve 
the vibrational displacements. Therefore, the light-element function < u s oi y «oo2z > m ay 
exhibit the large values compared with the other two. This indistinguishability gives rise to 
experimental ambiguity on this study somehow. 

Fig. 3 (a)-(c) represent continuous modulation of the function values. Portions of the 
curves with negative < u s pb x uooiy >, positive < u S B X uooiy >, and positive < u s oi y uoo2z > 
(or these with the reverse signs) correspond to the polar correlations as mentioned above. 
Therefore, lengths of the correlations can read from the figures. In (a)-(c), there are three, 
four, and three kinds of appreciable peaks of the functions, respectively. Their tail-to-tail 

— * o 

separations in sequence of the peak appearance from | R s |=0 are ~80, 30, and 30A for 
| R s |//[111], -80, 30, 50, 50A for | R s |//[110], and -80, 60, 60A for | R s |//[100] 
(< u sOly u o2z > of the second and third peaks are merged with each other. But because of 
a dip at the center, we realize that there are two different peaks. The other two peaks are 
separated clearly). These data suggest complicated three-dimensional spatial polar modu- 
lations. However, we may recognize from the sequence of the lengths that the correlations 
make spheres or ellipses whose diameters are 30-80A. These values agree with the length 
order found in the previous TEM [9] and neutron work [10] which reported that the size 
is ~100A and 30-40 A, respectively. Since TEM measures the length by image contrast, 
observed diffuse interfaces of the microregions may give this difference. On the other hand, 
the difference from the neutron result can be attributed to our technical indistinguishability 
of the dynamic displacements as mentioned above. 

Fig. 1 (c) shows the contour map given by the data taken at 270K, corresponding to 
the map at 300K. The results indicate increase of the intensities near the zone centers (The 
contours are significantly close to each other making dark areas.) with no essential difference 
in the distribution from that at 300K, which is also consistent with the previous x-ray work 

7 



[3]. Similarly, the observed data were analyzed quantitatively, and the recalculated map 
is represented in Fig. 1 (d). The observed data are well reproduced giving the R factor 
~3%. The three displacement-amplitude products, < u Pb £ x u 01 j: y >, < u B ^ x u Q1 ^ y >, and 
< u o\ky u 02kz > ' are a ^ so significantly larger than the others, and construct the same local 
symmetry as that at 300K. The correlation functions calculated are plotted in Fig. 3 (d)- 
(f), corresponding to (a)-(c). The results give significant reduction of main-peak values of 
the functions for the three ionic pairs by a factor of ~l/300, ~l/70, and ~1/1000. This 
drastic change suggests freezing of thermal lattice fluctuations [3] which may occur near 
temperature where a dielectric constant peaks, i.e., ~270K. The decrease of the thermal 
agitation allows new peak structures to appear in the function profiles giving approximately 
spherical correlations with a diameter of 30-40A which is most likely to be a size of the polar 
microregions. This value is consistent with that reported by the previous neutron work [10]. 

The data taken at 240K exhibit no significant difference in the intensity distribution from 
that at 270K. The analysis also gives no change in the local symmetry and its correlation 
size but shows slight reduction of the function values. This variation may simply be caused 
by q-dependent freezing [3] or decrease of energy of uncorrelated phonons, and is suggested 
to be two-phase behavior found by the NMR work [8]. In fact, the function peak structures 
are more obvious compared with those at 270K. 

Below 240K, enlargement of the microregions occurs [9]. Therefore, the diffuse intensities 
are expected to condense closer to the zone center. This increase must be remarkable in 
our unaccessible zone-center region, ~(0.01) 3 rX-u. 3 , and was obstruction to perform the 
measurement at the lower temperature. The higher instrumental resolution is required for 
the further investigation. 

The polar domains which appear through the local condensation exhibit the local polar- 
izations randomly along anyone of the <111> directions giving the macroscopic paraelectric 
state, and the polarizations must be in a sort of frustrating states under an electric field, 
which leads to the dielectric dispersion. A zero-frequency picture of the random polariza- 
tions gives an idea that the Coulombic field interactions among the microdomains have a 
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tendency to prevent their growth each other from the embryonic states. This effect requires 
undercooling to promote the further transition. The TEM study [9] and another previ- 
ous synchrotron x-ray work [f 9] determined that temperature of the initiation is 220-230K, 
meaning undercooling of several-ten kelvins. Local lattice distortion generated by the charge 
misfit and the size effect simply provides the preferential condensation sites heterogeneously, 
and causes a broad maximum of the dielectric constant. The dispersion is just a consequence 
of the electric-energetical growth barrier. 

N.T. thanks Y. Uesu, K. Fujishiro, N. Yamamoto, and G. Shirane for valuable private 
communications on this study. The acknowledgment is extended to K. Ohshima, K. Ohwada, 
and Y. Katsuki for their assistance in the measurement. 
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FIGURES 

FIG. 1. The Oth-(IOO) contour maps of (a) observed I(qi) at 300K and (b) its calculated result, 
and correspondingly (c) and (d) at 270K. The contour values are indicated in 10 4 e.u. for all pieces. 
The contour interval is 0.8 in the units. 

FIG. 2. Schematic representation of (a) < u pb ^ x u 01 ^ y >, (b) < u B ^ x u Q1 ^ y >, (c) 
^ u oikx u 02ky ^' an( ^ their resultant shifts given in terms of the symmetry Pm3m. Superpo- 
sition of all ionic shifts shows the correlation which is rhombohedral polar as illustrated in (d), and 
the symmetry corresponds to R3m. 

FIG. 3. The correlation functions of intercell distance R s at 300K along (a) the [111], (b) [110], 
and (c) [100] directions. Figures of the functions at 270K corresponding to (a)-(c) are given as 
(d)-(f). Values of the three functions, < u s pb x uooiy >, < u s B x uooi y >, and < u s oi y uoo2z >, are 
indicated in each figure by circles, squares, and triangles, respectively. 
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